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INTRODUCTION 


Phase-shift  (References  1  through  5)  and  ring-down  (References  4 
and  5)  methods  of  determining  the  reflectance  of  mirrors  in  a  high  Q 
optical  cavity  possess  high  sensitivity  and  precision.  In  maintaining 
this  precision,  the  errors  in  the  measurement  must  be  minimized.  The 
errors  include — but  are  not  limited  to — (1)  mode  matching  the  input 
source  to  the  optical  cavity,  (2)  stable  input  source  in  power  and 

waveform  over  the  time  required  to  make  measurement,  (3)  in  terms  of 
loss,  cavity  mirrors  equal  to  or  better  than  the  unknown  mirror's 

losses,  and  (4)  thermal  stability  of  the  cavity.  Rahn  (Reference  5) 
examined  an  optical  cavity  and  included  a  random  thermal  fluctuation 
in  the  analysis  to  explore  the  range  of  fluctuations  that  can  be 
tolerated.  He  was  not  able  to  evaluate  two  time-averaged  integrals 
exactly  but  derived  an  expression  for  the  loss  on  a  numerical 

approach.  Starting  from  the  same  assumptions  as  Rahn,  we  will  derive 
an  exact  solution  to  a  passive  optical  cavity,  including  the  thermal 
fluctuation.  Also,  an  error  analysis  of  the  phase-shift  method  is 

derived  and  examined. 


THEORY 


Consider  an  optical  cavity  consisting  of  two  mirrors  with  their 
symmetry  axes  on  a  common  line,  as  in  Figure  l.  Once  the  reflectance 
of  one  of  the  mirrors  is  known,  the  other  can  be  calculated.  A  phase 
shift  can  be  uniquely  associated  with  the  product  of  the  cavity  mirror 
reflectances: 


where  Rc  is  the  total  reflectance  of  the  cavity  and  Rj  and  R2  are 
the  respective  reflectance  values  of  the  cavity  mirrors. 


INJECTED 
LASER  SOURCE 


DETECTOR 


FIGURE  1.  Two-Mirror  Passive  Cavity.  Injection  laser  source  can  be 
short  pulsed  laser  or  modulated-CW  laser  source. 


Another  configuration — a  three-mirror  cavity — will  be  discussed 
in  the  remainder  of  this  paper.  In  a  three-mirror  passive  cavity,  the 
optimum  position  of  the  test  mirror  is  at  the  vertex  of  the  triangle 
defined  by  the  three  mirrors.  In  the  configuration  of  Figure  2,  the 
total  reflectance  of  the  cavity  is 


INJECTED 

LASER 

SOURCE 


MIRROR  1 


MIRROR  2 


MIRROR  3 


FIGURE  2.  Three-Mirror  Passive  Cavity.  Advan¬ 
tage  over  cavity  in  Figure  1  is  in  one-round 
trip  test  Mirror  3  is  sampled  twice,  i.e., 
Rt  =  R].R2^3^* 
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Rt  =  R1R2R3 


where  R3  is  the  test  mirror  reflectance.  R3  is  squared  because  in  one 
round  trip  the  test  mirror  is  sampled  twice.  Thus,  R3  may  be  calcu¬ 
lated  if  Ri,  R2,  and  Rt  are  known.  To  obtain  high-reflectance 
values  of  R3 ,  high-reflectance  end  mirrors  are  required  to  aid  in 
minimizing  the  error  of  the  measurement.  The  effects  of  errors  will 
be  discussed  later. 

To  model  the  cavity  exactly,  careful  attention  must  be  given  to 
the  input  wave  and  its  representation.  The  representation  of  the 
input  laser  power  used  by  Rahn  (Reference  5)  will  be  the  starting 
point  in  this  analysis.  This  power  is  described  by 


I(t)  =  I  sin2  wt  =  4  I  (1  -  cos  2<dt) 
o  2  o 


where  IQ  is  the  output  power  of  the  laser,  f  (a)  =  2  7tf)  is  the  ampli¬ 
tude  modulation  frequency  of  this  laser  determined  by  a  suitable 
intensity  modulator,  and  t  is  time. 

The  passive  cavity  is  designed  to  allow  oscillation  of  the  cavity 
modes  defined  by  fQ  -  nc/2L,  where  fQ  is  the  optical  frequency  of 
the  injection  source  (units  =  sec-1),  n  is  any  integer  (on  the  order 
of  10°  for  optical  frequencies  injected  into  cavities  with  lengths 
from  centimeters  to  meters),  L  is  the  length  of  the  cavity  in  meters, 
and  c  is  the  speed  of  light  in  meters  per  second.  If  the  source  is  a 
laser  possessing  a  mode  bandwidth,  a  change  in  cavity  length  will  not 
quench  the  oscillation  but  will  select  a  slightly  different  mode.  If 
the  length  of  the  cavity  changes  beyond  an  integral  number  of  half 
wavelengths  or  if  only  one  mode  is  available  from  the  laser,  the 
oscillation  may  be  intermittent.  This  will  lead  to  errors  in  measure¬ 
ment  of  the  cavity  loss.  Care  in  system  design  is  necessary  to  pre¬ 
vent  intermittent  operation. 

The  beam  is  free  to  propagate  through  the  cavity  subject  to 
losses  from  the  mirrors.  Since  the  modulation  of  the  beam  changes  its 
phase  with  time,  we  can  expect  that  as  the  beam  traverses  through 
the  cavity  in  the  course  of  one  round  trip,  the  relative  phase  changes 
by  the  amount  (References  6  through  8) 
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U)T  =  (2itf)( — )  = - 
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for  a  single  mode.  The  resulting  Lntercavity  intensity  is  therefore  a 
summation  of  the  Intensities  of  all  modes  oscillating  in  the  cavity. 
Each  of  these  modes  is  subject  to  reductions  in  intensity  because  of 
the  cavity  losses.  If  Equation  3  is  representative  of  the  modulation, 
this  summation  over  all  modes  can  be  expressed  as  (References  5  and  8) 


I  *  I  y  Rn  sin2  (ojt  +  nun) 
o  L  t 


where  i  **  2L/c  and  Rt  is  the  total  cavity  reflectance.  This 
summation  can  be  evaluated  by  expressing  the  trigonometric  function  in 
its  complex  exponential  form,  so  that  the  summation  takes  the  form  of 
a  geometric  series.  In  the  case  of  Equation  3,  the  appropriate 
trigonometric  power  relation  is  applied.  Details  of  this  process  are 
Included  in  Appendix  A.  Thus, 


00 

I  -  IQ  l  Rt  si °2  (wt  +  nwT)  "  1o  2U^rT 


(sin  2u>t)(R  sin  2u)t)-(cos  2<ot)(l-R  cos  2u>t) 
_ t _ t _ 

1  -  2R  cos  2wt  +  R^ 


The  power  and  the  phase  angle  of  the  transmitted  beam  It  can  be 
related  since  it  follows  a  pattern  similar  to  that  of  the  injected 
beam.  Thus , 


!T  a  -y  [!  "  cos  (2wt  +  <Ji)  ] 


4-  I  +  I  (sin  2 'tit  sin  <i>  -  cos  2 at  cos  <{>) 
2  o  2  o 


•v-v  «> 
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Now  the  time  dependence  can  be  eliminated  by  equating  the  appropriate 
time-related  coefficients  in  Equations  3b  and  4.  Thus,  equating  the 
coefficient  of  sin  2u>t, 


R  sin  2  arc 

sin  =  - 

1  -  2Rt  cos  2ut  + 


(5) 


Equating  the  coefficient  of  cos  2<i)t  in  Equations  3a  and  4, 


l  -  R  cos  2oji: 

cos  <t>  - - 

1  -  2Rt  cos  2ojt  +  R^ 


(6) 


So  that 


tan 


sin  4) 

COS  (j) 


R^  sin  2(i)T 
1  -  Rt  cos  2ut 


(7) 


This  is  the  result  if  no  random  thermal  fluctuations  are  included  in 
the  analysis. 

To  examine  the  effect  of  a  random  phase  induced  by  a  thermal 
variation  in  the  cavity,  the  amplitude  rather  than  the  intensity  must 
be  used.  The  amplitude  function  is  related  to  intensity  I  by 


I  -  A*A  (8a) 

as  is  also 

Rt  =  r*r  (8b) 


where  A  is  the  electric  field  amplitude  of  the  laser  beam,  r  is  the 
reflectivity  coefficient,  and  Rt  as  before  is  the  reflectance  of  the 
cavity. 

Rahn  (Reference  5)  suggests  that  a  random  phase  shift  6  could 
range  from  -n  to  +x;  he  used  numerical  averaging  to  obtain  an  expres¬ 
sion  for  <tan  $>.  Increased  precision  may  be  realized  if  the 
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averaging  evaluation  of  tan  <}>  can  be  done  exactly.  Starting  with 
Equation  4  of  Reference  5  for  the  transmitted  intensity. 


A(t)  =  A(jt  jt  2t  3  (r)F  sin[io(t  -  p  x )  J 

p=0 


where  r  -  rir2r3  »  tj,  t2,  and  tj  are  the  reflection  and  transmission 
coefficients,  respectively;  p  Is  the  number  of  passes  through  the 
cavity;  6  is  the  round-trip  phase  shift;  and  x  is  the  time  required 
for  a  photon  to  make  a  single  round  trip  in  the  cavity.  Summing  the 
Infinite  series  gives 


sln  _ \ _ . _ i _ 

M  2  2  ,  i(8  P  i.tx)  ,  1(8  -  u>x) 

l  -  re  1  -  re 


1  cos  bit _ 1 _ _  1 _ 

2  i(5  +■  >ox)  .  i(6  -  'ox) 

1  -  re  1  -  re 


sin  wt(l  -  r^je^  cos  <ox)  -  cos  wt(rir2e1^  sin  ux) 


l  -  2rel°  cos  <ox  +  r  e 


202i  6 


Squaring  to  get  the  transmitted  power. 


IT  -  | A(t) 


1  1  2 

cos  2uit(r  cos  iox  cos  6  -  x  -  -=  r  cos  2oox) 

2  _  _  2  L _ 

2  2  /  4 

4r  cos  b  -  4  cos  <ot(  r  +  r  3)  cos  6  +  4r 


-  sin  2(ot(r  sin  t ox  co s  <S  -  sin  urt  cos  (ox) 

-  2r2  +  4r2  cos 2  a>x  f  1 


details  of  which  are  in  Appendix  B.  Equating  trigonometric 
coefficients  in  Equations  4  and  13  for  the  transmuted  power  gives 
expressions  for  sin  d>  and  cos  0, 


s& 
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2 

2(r  sin  cot  cos  cox-r  sin  cox  cos  6) 


4r  cos  6-4  cos  cox(r+r  )  cos  6+4r  -2r  +4r  cos  cox+1 


(1  +  r^  cos  2cox)  -  2r  cos  cox  cos  6 


4r  cos  6-4  cos  oox(r+r  )  cos  6+4r  -2r  +4r  cos  cox+1 


Both  expressions  are  of  the  form 


a  cos  6  +  b 


c  cos  6  -  d  cos  6  +  e 


where 


c  =  4r‘ 


d  =  4r(l  +  r  )  cos  cox 

5  o  4  0 

e  »  1  +  4r  cos  cox  +  r  -  2r 


e  =  (1  -  r^)2  +  4r2  cos 2  cox 


For  sin  <t>. 


a  =  -2r  sin  cox 


b  =  2r  sin  cox  cos  cox 


For  cos  4>, 


a  =  -2r  cos  cox 


b  =  1  +  r  cos  2  cox 


Equations  11  and  12  rewritten  in  the  form  of  Equation  13  must  be 
integrated  between  the  limits  0  to  2u  to  obtain  an  average  value  of 
tan  <fc.  Evaluation  of  these  integrands  is  by  partial  fractions.  The 
decomposition  of  the  denominator  to  produce  partial  fractions  is 
accomplished  by  noting  that 
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c  cos 1  6  -  d  cos  6  +  e  =  (/c  cos  6  -  h)(/c  cos  6  -  k) 

r\  - - 

=  c  cos  6  -  (h  +  k)/c  cos  6  +  hk 

where  h  and  k  are  coefficients  that  satisfy  d  =  (h  +  k) /c  and  c  =*  hk. 
These  equations  yield 


and 


k/c  -  kd  +  e/c  =  0 


-  k 


so  that 


1 


1 


k  =  — - —  Id  +  /d c  -  4ec  )  ,  h  =  — (d  -  /d  ^  -  4ec 

2/c  \  /  2/c  \ 


To  find  the  partial  fractions,  we  first  write 


a  cos  6  +  b 


+ - 


B 


c  cos^  6  -  (h+k)  /c  cos  b  +-  hk  /c  cos  6  -  h  /c  cos  6  -  k 


where  A  and  B  are  unknown  coefficients,  then 


a  cos  6  +  b  =  A(/c  cos  <S  -  k)  +  B( /c  cos  6  -  h) 


and  evaluate  at  /c  cos  6  =  k  and  /o  cos  6  =  h  separately.  These  yield 
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-(ah  +  b/c) 


( 14a) 


ak  +  b/c 
/d2  -  4ec 


(14b) 


We  integrate 


.  2rc  (a i  cos  9  +  bj)d 9 

<sin  <t>>  =  -^—r  J  - 2 - 

Z71  o  c  cos  9  -  d  cos  9  +  e 


)n  d9 

2 11  o  /c  cos  9  -  h 


,  Bl  2fU  d 9 

+  J  _ 


Al  BX 

-  +  - 


o  /c  cos  9  -  k 


A2  - 


(15a) 


with  similar  results  for 


^  2n  (a2  cos  9  +  b2)d9 

<cos  *>  =  ^  J  - 2 - 

o  c  cos  9  -  d  cos  9  +  e 


(15b) 


so  that 


<tan  *>  ,  <i±ij*>  , 


Ai/k2  -  c  +  Bi/h2  - 


<C°8  ♦>  A2A2  -  c  ♦  B2X2  - 


This  becomes 


r  sin  2u)X 


<tan  <{>>  = 


Rt  sin  2  art 


,  2  0  1  -  R_  cos  2wt 

1  -  r  cos  2wx  t 


1 1 


**•1 
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the  details  of  which  appear  in  Appendixes  C  and  D.  Notice  this  is 
exactly  equivalent  to  Equation  7.  Thus,  random  phase  shifts  are  not  a 
major  source  of  error  over  the  time  interval  of  measurement  of  <J>. 

An  expression  for  R3 — the  unknown  mirror's  reflectance — can  be 
derived  starting  from  Equation  17.  Using  a  common  denominator  and 
solving  for  Rt  gives 


t  _  _ <tan  4>> _ 

t  cos  2tox  <tan  <)>>  +  sin  2 «x 


Using  Equation  2,  solving  for  R3,  and  substituting  into  Equation  18 
for  Rt  gives 


_ <tan  <ft> _ 

RjR2(cos  2cjx  <tan  $>  +  sin  2wx) 


The  cavity  parameters  4xfL/c  substitute  for  wx;  thus 


_ <tan  <b> _ ' 

n  „  8xfL  .  ,  8xf L 

r1r2  Ktan  d>>  cos  — ^ —  +  sin  — - — 


If  the  cavity  satisfies  8xfL/c  <1  and  R^  ~  R2  ~  1,  then 


/  <tan  (j)> 

3  is*.  ^  ,  8xfL 

\  <tan  4>>  +  — - — 


a  good  approximation  for  very  low-loss  mirrors. 


A  comparison  of  the  sample  reflectance  R3  derived  by  the 
numerical  analysis  of  Rahn  (Reference  5)  and  the  exact  expression, 
Equation  20,  is  given  in  Figure  3.  The  expression  used  by  Rahn  is 
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<tan  0>  -  4xfL 


RiRo  <tan  0> 


EXACT // 

// 

yx  NUMERICAL 


It® 

IW 
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& 


PHASE  ANGLE  (J>,  DEG 

FIGURE  3.  Comparison  of  the  Numerical  Analysis  by  Rahn  (Reference  3) 
and  the  Exact  Analysis  for  the  Same  Cavity  Parameters,  Ri  *  R2  = 
0.998,  f  =  100  kHz,  and  L  =  3.864  ra.  Shows  sample  reflectance  versus 
measured  phase  angle  between  injected  and  transmitted  wave. 

where  L  is  the  round-trip  length  of  the  cavity.  Recall  that  the 
length  used  in  Equation  20  is  defined  as  the  length  of  the  cavity. 
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Notice  that  the  numerical  analysis  gives  a  lower  sample  reflectance 
than  the  exact  expression,  but  differences  are  not  seen  in  R3  until 
the  fourth  decimal  place.  The  numerical  analysis  of  Rahn  is  adequate 
unless  one  requires  higher  precision,  1  part  in  104  or  better.  This 
is  the  situation  for  characterization  of  state-of-the-art  low-loss 
mirrors  used  in  laser  gyros  and  low-gain  lasers.  The  high  precision 
is  also  useful  for  detecting  very  small  changes  in  the  performance  of 
an  optical  coating  under  laser  loading  when  determining  functional 
lifetime. 


ERROR  ANALYSIS  FOR  THREE-MIRROR  CONFIGURATION 


In  order  to  evaluate  the  total  error  in  the  value  of  reflectance 
of  a  sample  mirror  in  a  passive  cavity,  errors  must  be  included  from 
four  sources:  errors  from  (1)  modulation  frequency;  (2)  cavity 
length;  (3)  cavity  end-mirror;  and  (4)  phase  in  the  phase  detector. 
Although  several  of  these  contribute  relatively  insignificantly  to  the 
overall  error,  the  error  analysis  nevertheless  includes  them,  both  as 
reference  for  future  systems  and  for  the  sake  of  completeness.  Future 
advances  in  measuring  systems  may  reduce  some  errors  to  the  point  that 
others  are  increasingly  significant.  Existing  systems  can  have  their 
total  errors  calculated  without  concern  that  error  sources  have  become 
more  (or  less)  significant. 


tan  4>  = 


Rt  sin  (2wt) 

1  -  R^  cos  (2  ore) 


So  that 


a*  <tan  = 


sin  (2uit) 

1  -  R  cos  (2gjt) 


( 1— R  cos  2orr)  sin  2u>x+R_  sin  2ort  cos  2orc  dR 
2,.>  t  t  t 

sec  ($)  =  - - - -t— 

(1  -  R  cos  2o>t)  2  0<5> 
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R  cos  2u>x 


(1  -  Rt  cos  2ujx)  sin  2u>x  +  Rfc  sin  2u>x  cos  2u>x 


1  -  R^  cos  2wx 


esc  (2oix) 


Since  we  are  concerned  with  the  error  associated  with  the  sample 
mirror  rather  than  that  of  the  system,  we  must  consider  the  previous 
equation  Rt  =  RlR2R32»  where  Rj  and  R2  are  the  system  mirrors  and  R3 
is  the  sample,  so  that 


✓RLR2  2/5^ 


1  t 

2R1R2R3  9$ 


The  chain  rule  applies  for  the  variables  in  the  equation  relating 
to  the  frequency  and  delay  time  of  the  system;  also,  so  that  the  total 
error  is 
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Evaluating  the  partial  differentials  gives  the  following: 


9  _ tan  4> _ 

3x  sin  2o)t  +  cos  2fjx  tan 


-  tan  <ft(2u)  cos  2u)x  -  2co  sin  2o)X  tan  d>) 

2 

(sin  2ojx  +  cos  2ojt  tan  4>) 

-2oi  tan  4>(cos  2<ox  -  sin  2u)x  tan  <t>) 

2 

(sin  2wx  +  cos  2ojx  tan  <{ i) 


_ tan  <6 _ 

sin  2wx  +  cos  2wx  tan 


-  tan  4>(2t  cos  2u)t  -  2  x  sin  2o)x  tan  <j>) 

2 

(sin  2u)x  +  cos  2wx  tan  <j>) 

-2x  tan  $  (cos  2u>t  -  sin  2 cox  tan  $) 

2 

(sin  2ojx  +  cos  2wx  tan  4>) 


2  2 
A  tan  $(cos  2ux  -  sin  2u>x  tan  <J>)  r/  .  -2.,  »  \2, 

- — - [(x»Au)  +((!>•  Ax)  ] 

(sin  2u>x  +  cos  2u>x  tan  $)  ** 


i  J  M  S)  «  m'<  j'*  * '  »  #  Ii*  ft*1  Jt1"  *  ft  >4 '  ul'kl'ulVI'i'  I  **  *  j  * .  *.» 
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4Rt2  sin2  2ort 


cos  2uit  - 


Rt  sin2  2  art  \2 


1  -  Rt  cos  2ort 


(1  -  R  cos  2  art)' 


(R  sin  2  art  cos  2 ojxY4 

sin  2urt  +  - 5 - - - ^ - ) 

^  1  -  R^cos  2ort  J 


*  [(  x*Aoi)  2  +  (oj*At)  2J 


4R2  sin2  2ort(l  -  R  cos  2ort)2 


(1  -  R  cos  2  art) 


cos  2  a rt  - 


2 

R  sin  2  art 


1  -  R  cos  2  art 


sin  2art  + 


R  sin  2urr  cos  2  art 


i  -  Rt  cos  2art 


[(x»Aw)2  +  (a>*Ax)2] 


4R2[sin  2wx  cos  2art(l-Rt  cos  2urt)-Rt  sin3  2art]  2[(  x»Aoj)  2+(  a>*Ax)  2] 


[sin  2ojx(  1-R  cos  2  art)  +  R  sin  2ojx  cos  2  art] 


4R2  [sin  2 art  cos  2  art  -  R^  sin  2u>x]  2  [(x»Aa>)2  +  (a>»Ax)2] 


sin  2arr 


4R2  (cos  2a)x  -  Rj.)2  t(x*Aa))2  +  (a)»Ax2] 


sin  2a >x 


Aw  =  A(2itf)  -  2 nAf  ;  Ax  =  A(— )  *  —  AL 
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2  2 
sec  4>  =  1  +  tan 


so  that  Equation  22  becomes 


1  —  2  cos  2ujt  +  R* 
(l  -  R£  cos  2wt)2 


dR^  (1  -  R  cos  2ut)2  (1  -  2R^  cos  2u>t  +  Rt)2 


^  sin2  2ojt(1  -  R  cos  2  art) 2 


(I  -  2R^  cos  2o)x  +  R^)' 
2 

sin  2gjt 


According  to  Equation  2,  the  reflectance  of  the  sample  mirror  can  be 
expressed  as 


where  Rt  is  the  reflectance  of  the  three-mirror  system  and  Rj  and  R2 
are  the  end-mirror  reflectances.  Since  they  are  inextricably  bound  to 
the  measuring  system,  the  end  mirrors  contribute  to  errors  in  the 
determination  of  the  sample  reflectance  as  follows: 


dR3  /R  -  _  1/  2  /R„ 

nr-— 4-  i  =  — 

dRl  ^  6Rl  /Ri 


1  *\ 


/R^R2 


(27a) 


(27b) 
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so  that  the  total  squared  error  contribution  from  end-mirror  accuracy 
is 


il“>)  +(SiM2) 


,  2  2  2  2, 

x  [R2AR|  +  R1AR2] 


2  /  2  2 
1 


4  3  V  2  2  /  222 

\  R1  r2/  4R !R2R3 


and  the  total  root-mean-square  (rms)  error  is 


AR3  *  2(R1R2R3)  ~  2Rt  COS  2“T  +  r2)csc2(2wt)A4>] 


+  4R2(cos  2arr  -  Rt)  2  esc2  iot 


x  4[(xAtf)  2  +  (-^)2J  +  R34(R22AR12  +  R12AR22)  (29) 


For  example,  a  system  In  which 


f  -  105  ±  102  hertz 


L  *  2  meters  ±  0.002  meters 


Rl  -  R2  -  0.999  ±  0.0001 


then  o)T  equals 


art  =  — ■  0.48  degrees 


A  sample  with  a  nominal  reflectance  of  0.998  (therefore,  Rt  *  0.996) 
would  produce  an  angle  exiting  the  cavity  of 


R  sin  2 an; 

$  -  tan-1  — - - - r —  =  82.706  degrees 

Y  1  -  Rfc  cos  2ojt  0 
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If  this  measurement  is  subject  to  an  error  of  +1  degree,  then  the 
reflectance  measurement  is  subject  to  an  error  of 


AM*)  “  2r7^r7  (1  ~  2Rt  cos  2a,T  +  Rt2)  csc  (2wT)(lfe) 


=  1.545  *  10- 


The  measurement  due  to  end-mirror  deviations  would  be 


AR3(  R  iR2) 


2R1R2R3 


[T  222  r3  f~i  2 

/R2AR1  +  R1AR2  =  2R ^~R 2  ,'2R1^R1 


R3ARi 
✓2  Ri 


7.064  x  10- 


Finally,  the  error  caused  by  the  frequency  and  cavity-length 
deviations  would  be  equal  to 


AR3(  wx) 


'( icf  Ax) 


2  .  ,LAu) 


R1R2R3 


*<■*?>  (  \  \  . 

2R3  U  -  Rt  cos  2wx J 


cos  2wx  -  Rt> 


/(4.1888xl0-V.(4.1888x,l0-t)2  (240.9)(3.85xI0-3) 


=  5.516  x  10* 


80  that  the  total  rras  error  would  be  equal  to 
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AR3  fc  =  /(l . 545  x  10"4)2  +  (7.064  x  10~5)2  +  (5.516  x  10"6)2 
tot 


1.700  x  10' 


and  the  reflectance  of  the  sample  would  be  expressed  as  R3  =  0.998  ± 
0.00017.  It  is  seen  from  the  above  that,  for  this  system,  the  major 
sources  of  error  in  calculating  the  reflectance  are  the  measurement  of 
the  phase  angle  of  the  superimposed  waves  and  the  uncertainty  of  end- 
mirror  reflectances.  However,  the  contributing  sources  of  error  are 
sufficiently  comparable  in  magnitude  that  we  suggest  including  them 
all,  especially  with  the  foreseeable  modifications  of  greater  phase- 
angle  sensitivity  and  different  end-mirror  uncertainty  values.  The 
above  analysis  would  not  be  valid  for  a  two-mirror  system,  since 
Rt  ■  RiR2»  but  those  elements  that  are  based  on  partial  derivatives 
of  Rt  can  be  incorporated  into  a  two-mirror  system  error  analysis. 


SUMMARY 


An  exact  calculation  has  been  performed  on  the  influence  of  a 
random  phase  fluctuation  induced  by  a  thermal  variation  in  the  cavity 
to  the  performance  of  a  passive  cavity  ref lectometer.  The  analysis 
started  with  the  assumptions  stated  by  Rahn  (Reference  5)  in  his 
numerical  analysis  of  this  problem.  The  expression  that  was  derived 
is  identical  to  the  expression  where  no  random  phase  fluctuations  were 
included.  Thus,  thermal  variation  of  the  cavity  is  not  a  major  source 
of  error  for  the  operation  of  a  ref  lectometer.  Comparison  of  the 
numerical  and  exact  analyses  show  that  the  numerical  expression  for  R3 
is  not  a  good  approximation  where  high  precision — better  than  1  part 
in  104 — is  required. 

Error  analysis  was  performed  to  determine  the  functional 
dependence  of  cavity  length,  modulation  frequency,  cavity  mirror 
reflectance,  and  phase  measurement  errors  on  the  precision  of  the 
ref lectometer.  The  total  error  of  a  simple  ref lectometer  was 
determined  to  be  on  the  order  of  2  parts  in  10 4 . 

If  high  precision  is  required,  care  must  be  taken  in  the  design 
and  construction  of  a  passive  cavity  ref lectometer  to  minimize  or 
eliminate  the  sources  of  error  that  are  mentioned  above.  On  the  other 
hand,  a  simple  ref lectometer  can  be  constructed  with  a  precision  and 
error  of  approximately  2  parts  in  104;  either  analysis  is  adequate  for 
Che  measurement  of  R3. 
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Appendix  A 

DERIVATION  OF  THE  INFINITE  SUM  IN  EQUATION  3 


S 

s 

ft 


Using  the  power  relation  for  sinz  (a), 


OO  00 

I'  «*  I  sl°2  +  nwt)  =  J  £  R^[l  -  cos  (2wt  +  2nurt)] 


Now  applying  the  sum-and-dif ference  formula 


OO  00 

I'  »  -i  £  Rt  ~  Y  ^  R"[cos  2wt  cos  2nwt-sin  2wt  sin  2nurc]  (A-l) 


I'  =  j  (y~  j  ^-)  -  cos  2a)t  l  R^  cos  2nu>r 


OO 

+  -j  sin  2wt  £  sin  2ni 


(A-2) 


The  trigonometric  identities  for  complex  exponential  yield 
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co  R 

Tl  1  1  „  v  t  ,  i2nojt  ,  -i2nu>TN 

1  =  in  -~ry "  i cos  2,ot  i  t  (e  +  e  > 

t  n=o 


+  ~  sin  2o>t  \ 


®  i2nujx  -i2no)x 

r  „n  e  -  e 


1  1  „  „  r  „n  i2niox  r  Dn 

t  r~R~)  "  -4  cos  2a)t  l  Rte  +  l  Rte 


n  -12no)X 


20  -  Rt) 


1  r  „n  i2nujx  r  „n 

+  —  sin  2wt  )  R  e  )  R„e 

4i  '■  t  L  t 


n  -i2ncox 


2(1  -  Rt)  4 


COS  2<dt 


- ! - + _ ! - 

1  -  R^12”1  1  -  Rte-12“T 


+  —  sin  2wt  - T7. - T7i - 

41  .  D  i2u>x  _  -12u)x 

1-Re  1-Re 

t  t 


(A-3) 


The  complex  fractions  can  be  evaluated  by  expressing  them  with  a 
common  denominator: 
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Appendix  B 


A*A 


DERIVATION  OF  A*A 


Given 


A  ~  sin  cot 


1  -  re  cos  cox 


.  „  16  2  2i 6 

I  -  2re  cos  tox  +  r  e 


+  cos  (Ot 


16  , 

re  sin  cox 


,  «  i6  2  216 

1  -  2re  cos  cox  +  r  e 


i  6  i  6 

sin  (Pt(  1  -  re  cos  cox)  +  cos  (ot(re  °  sin  cox) 

,  o  15  ,  2  2i6 

I  -  2re  cos  cox  +  r  e 


sin2  cot  (1-re1^  cos  cox)(I-re  ^  cos  tox)+cos2  cot(r2  sin 


, ,  .  i6  2  2i6,,  ,  „  -16  2-2 

(1  -  2re  cos  cox+re  )(1  -  2re  cos  cox  +  r  e 


+  sincotcoscot  ( 1-re1  ^coscox)re  1 ^sincox+(  1-re  1  ^coscox)re1  ^s 


-i  6 


-i  6 


i6 


/.  0  15  2  216.,.  -i6  2  -2i6  . 

(l-2re  coscox+r  e  )(l-2re  coscox+r  e  ) 


The  denominator  becomes 
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.  -  -i6  ,  2  -2i6  „  i6 

1  -  2re  cos  arc  +  r  e  -  2re  cos  wx 


,2  2  „  316 

+  4r  cos  wx  -  2r  e  cos  wx 


2  2i6  0  3  i6  4 

+  r  e  -  2r  e  cos  wx  +  r 


4  2  2  2 

1  +  r  -  4r  cos  6  cos  cjt  +  2r  cos  26  +  4r  cos  ojt 


-  4r  cos  6  cos  wt 


4  2  2  2 

=1  +  r  -  4r(l  +  r  )  cos  6  cos  wx  +  4r  cos  wx 


+  2r 2(2  cos2  6-1) 


2  2  2  2  2 
=  4r  cos  6  -  4r(l  +  r  )  cos  6  cos  wx  +  4r  cos  wx 


+  (1  -r2)2 


Simplifying  the  numerator  as  follows. 


|  (1  -  cos  2o>t)(l  -  re  ^  cos  ojx  -  re  1  ^cos  wx  +  r2  cos2  wx) 


+  ^  (1  +  cos  2wt)(r2  sin2  wx) 


+  sin  2wt(re  ^  sin  wx  -  r2  sin  wx  cos  wx 


i6  .  2  .  x 

+  re  sin  wx  -  r  sin  wx  cos  wx) 


1  2  2  2  2 
Y  ( 1  -  2r  cos  6  cos  wx  +  r  cos  wx  +  r  sin  wx) 


It  is  noted  that  the  expression 


/<1  ^  —  4ec 

is  an  imaginary  term  since 

d2-4ec  *  16r2(l+r2)2  cos  2  wt-4[(  1-r2)  2+4r2  cos  2  u)x](4r2) 
*  -16r2(l-r2)2  sin2  u>x 


Thus  -  4ec  "  4ir(l  -  r2)  sin  wx,  and  imaginary  terms  are  present 
in  Equation  16.  Upon  multiplying  the  numerator  and  denominator  of 
Equation  16  by  47"-  c. 


Ax/(k2  -  c)(h2  -  c)  +  BL(h2  -  c) 

<tan  4»>  =  - -  --  . . . —  - -  (C-2) 

A2/(k5  -  c ) ( h 2  -  c)  +  B2(h2  -  c) 

2  2  2 

Using  hk  *  e  and  (h  +  k  )c  *  d  -  2ec,  the  argument  under  the  radical 
in  Equation  C-2  becomes 
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(k2-c)(h2-c)  -  h2k2-(h2+k2)c+c2 


2  9  2 

*  e  -d  +2ec+c 


(e+c) 2-d2 


[( 1-r  2)  2+4r  2co8  2ojx+4r  2]  2-16r  2(  1+r  2)  2cos  2wx 


[(l+r2)2-4r2  cos2  wx]  2 


(C-3) 


X  *  /(k2  -  c)(h2  -c)*(l+r2)2  -  4r2  cos2  wx  is  entirely 
Now  Equation  C-2  can  be  written 


AXX  +  -  c) 

<tan  4 1>  »  - - -  ,  X  *  (1  +  r2)2  -  4r2  cos2  wx 


A2X  +  B2(h^  -  c) 


a  jh  +  b^/c  a  jk  +  bj^/c 

-  X  +  -  (h2  -  c) 

(h  -  k)*^c  (k  -  h)/c 

r-"  1  — -  1  'i  -  ■  ■  -  -  —  ■  ■  -  —  ■  ■  ■■  —  —  ■ 

'  a  2h  +  b  a  2k  +  b  2/c 

-  X  +  - 3-  (h2  -  c) 

(h  -  k)/c  (k  -  h)/c 


(ajh  +  b^/c)  X  -  (ajk  +  b^/c)(h  -  c) 
(a2h  +  b2/c)  X  -  (a2k  +  b2/c)(h2  -  c) 


(ajh  +  bj/c)  X  -  (ajkh2  -  a  jkc  +  bj/ch2  -  bjc/c) 
(a2h  +  b2/c)  X  -  (a2kh2  -  a2kc  +  b2/ch2  -  b2c/c) 


(C-4) 
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Using  hk  *  e  again. 


<tan  <p> 


(aih  +  bj/c)  X  -  (ajeh  -  ajkc  +  b^/ch2  -  b^c^c) 
(a2h  +  b2»^c)  X  -  (a2eh  -  a2kc  +  b2*^ch2  -  b2C»^c) 


ajh(X  -  e)  +  b]/c(X  +  c)  -  b^/ch  +  a^kc 
a2h(X  -  e)  +  b2*^c(X  +  c)  -  b2/ch2  +  a2kc 


ajh(e  -  X)  -  bj/c(X  +  c)  +  b^h  /c  -  ajkc 
a2h(e  -  X)  -  b2»^c(X  +  c)  +  b2h2/c  -  a2kc 


Using 


—  ( d  -  /d 2  4ec  ) 

!i^c  \  / 


— —  (d  -  Dl)  ,  D  ■  4r(l  -  r2)  sin  wx 
2/c 


h2  -  -p-  (d2  -  2dDi  -  D2) 
4c 


k  -  -±-  (d  +  Dl) 
2^c 


and  Equation  C-5  can  be  written 


(C-5) 
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al  _  bl  a  i  /c 

-  (d-Di)(e-X)-b]/c(c+X)+  — =  (d2-2idD-D2) - (d+Di) 

v  2fc  4/c  2  ,  „ 

<tan  $>  »  - — - (C-6) 

a  2  b2  a2»|,c 

-  (d-Di)(e-X)-b2Vc(c+X)+  -  (d2-2idD-D2) - (d+Di) 

2/c  4  7c  2 


or,  multiplying  numerator  and  denominator  by  4/c, 


<tan  <J>> 


2ai(d  -  Di)(e  -  X)  -  4b!c(c  +  X) 
2a2(d  -  Di)(e  -  X)  -  4b2c(c  +  X) 


+  bi(d2  -  2idD  -  D2)  -  2a!C(d  +  Di) 
+  b2(d2  -  2idD  -  D2)  -  2a2C(d  +  Di) 


(C-7) 


Separating  into  real  and  imaginary  parts. 


2aid(e-X)  -  4bjc(c  +  X)  +  bx(d2  -  D2)  -  2alCd 


<tan  <J)> 


2a2d(e-X)  -  4b2c(c  +  X)  +  b2(d2  -  D2)  -  2a2cd 


Di[2a^(e— X)  +  2b ]d  +  2a ^c] 
Di[2a2(e-X)  +  2b  2<i  +  2a2c] 


(C-8) 


A  =  2ajd(e  -  X)  -  Abjc(c  +  X)  +  bj(d2  -  D2)  -  2ajcd 
B  «  4r(l  -  r2)  sin  on 

C  =*  2a2d(e  -X)  -  4b2c(c  +  X)  +  b2(d2  -  D2)  -  2a2cd 

E  ■  2a2(e  -  X)  +  2b2d  +  2a2c 

G  =*  2a  i(e  *•  X)  +  2b  ^d  +  2a 
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<tan  $> 


A  -  BiG 
C  -  BiE 


.[AjL_BiG]  x  [C_L 

[C  -  BiE]  [C  + 


AC  +  B2EG  +  Bi(AE  -  CG) 


C2  +  B2E2 


(C-9) 


Since  the  original  integrand  was  real,  we  expect  that  the  above 
expression  should  also  be  real,  with  the  inference  that  AE  -  CG  *  0  or 


(C-10) 


This  is  indeed  the  case,  as  shown  in  Appendix  D.  Using  Equation  D-l, 


<tan  d>> 


A(^)  +  B2EG 

\f  )2  +  b2e2 


A2EG  +  B2EG3  =  EG(A2  +  B2G2)  _  G 
A2E2  +  B2E2G2  E2(A2  +  B2G2)  e 


(C-ll) 


ES 


m 

:: 

1 

k 

I 

I 

I 

I 

l 

i 

I 


k 


& 


J*  .li'.i* 


«*»','»  .•«  |*a  ■***  #»*'  **ir«*4<  Ji'.fk 


I »••  ***  »•!  a 


NWC  TP  6776 


<tan  $> 


2ai(e  -  X)  +  2b  id  +  2a ic 
2a2(e  ~  X)  +  2b  2d  +  2a 2c 


ai(e  -  X  +  c)  +  bid 
a2(e  ~  X  +  c)  +  b2d 


2  2  2  2 
2rsinu>x(  8r  cos  oix)-2r  6lmoxcosu)x(4rcos  wx)  ( 1+r) 


-2rcoswx(8r 2cos  2u>x)+(  1+r  2cos2t»)x)(4rcoswx)(  1+r2) 


3  2  3  2  2 

16r  sin  ojt  cos  ojt  -  8r  sin  wx  cos  wx(l  +  r  ) 


2  2  3  3 

4r  cos  wt(1  +  r  )(l  +  r  cos  2u>x)  -  16r  cos  wx 


3  2  2 

8r  sin  oix  cos  u)x(l  -  r  ) 


2  2  3  3 

4r  cos  wx(l  +  r  )(l  +  r  cos  2wx)  -  16r  cos  cox 


(C-12) 
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Dividing  the  numerator  and  denominator  by  4r  cos  cut. 


<tan  $> 


2  2 
_ 2r  sin  a>x  cos  arc(  1  -  r  ) _ 

(1  +  r2)(l  +  r2  cos  2  art)  -  4r2  cos2  art 


2  2 

_ 2r  sin  ux  cos  arx(  1  -  r  ) _ 

2  2  4  2  2 

1  +  r  cos  2u>x  +  r  +  r  cos  2<ox  -  4r  cos  u>x 


2  2 

_ 2r  sin  arc  cos  oit  (1  -  r  ) _ 

2  2  2  4  2  2 

1  +  r  (2cos  art  -  1)  +  r  +  r  cos  2arc  -  4r  cos  a>x 


<tan  <» 


<tan  $> 


2  2 
2r  sin  oit  cos  <ox(  1  -  r  ) _ 

1  -  2r2  cos2  art  +  r4(2cos2  arc  -  1) 


2  2 
2r  sin  oit  cos  art(  1  -  r  ) _ 

1  -  2r2  cos2  art  +  2r4  cos2  a>t  -  r4 


2  2 
2r  sin  arc  cos  a)t(l  -  r  ) 

(1  -  r4)  -  2r2  cos2  urt(l  -  r2) 


2  2 
_ 2r  sin  ait  cos  arc(  1  -  r  ) _ 

(1  -  r2)(l  +  r2)  -  (1  -  r2)(2r2  cos2  art) 


2r  sin  art  cos  an 
1  +  r2  -  2r2  cos2  arc 


r  sin  2  an  „  _  2 

- - -  ,  R  =  r 

1  -  r  cos  2ajx 


Rt  sin  2 art 
1  -  Rt  cos  2 art 


(C-13) 


(C-14) 
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Appendix  D 


PROOF  OF  AE  -  GC  -  0 
Starting  with  substitution  of  terms 
-  GC  -  [2aid(e  -  Xl)  -  Ab^  +  X)  +  b^d2  -  B2)  -  2aicd] 

x  [2a22a2(e  -  X)  +  2b2<l  +  2a2<:]  -  [2ai(e  -  X)  +  2b^d  +  2aic) 
x  [2a2d(e  -  X)  -  4b2c(c  +  X)  +  b2(d2  -  B2)-  2a2cd] 

-  4aia2(e  -  X)  2d  +  Aa^dCe  -  X)(b2d  +  a2c)  -  8a2  bjc(c  +  X) 
x  (e  -  X)  -  8bic(c  +  X)(b2d  +  a2c)  +  2a2bi(e  -X)(d2  -  B2) 

+  2bi)d2  -  B2) (b 2d  +  a2c)  -  Aa^cdCe  -  X)  -  Aa^d 
x  (b 2d  +  a2c)  -  (4aja2(e  -  X)  2d  -  8aib2c(c  +  X)  +  2aib2 
x  (e  -  X)(d2  -  B2)  -  4aja2cd(e  -  X)  +  4a2bid(e  -  X) 

-  8bib2cd(c  +  X)  +  2bib2d(d2  -  B2)  -  4a2bicd2  +  Aa^cd 
x  (e  -  X)  -  8aib2c2(c  +  X)  +  28^20(8*  -  B2)  -  4aia2c2d] 

■  4aid(e  -  X)(b2d  +  a2c)  -  8a2bic(c  +  X)(e  -  X)  -  8bic 


x  (c  +  X)(b2d  +  a2c)  +  2a2bi)e  -  X)(d2  -  B2)  +  2bx(d2  -  B2) 
x  (b2d  +  a2c)  -  4aia2cd(e  -X)  -  Aa^cd(b2d  +  a2c)  +  8aib2c 
x  (c  +  X)(e  -  X)  -  2a ib2(e  -  X)(d2  -  B2)  +  4a!a2cd(e  -  X) 

-  4a2bid2(e  -  X)  +  8bib2cd(c  +  X)  -  2bib2d(d2  -  B2) 

+  4a2blcd2  +  8aib2c2(c  +  X)  -  2aib2c(d2  -  B2)  +  4a1a2c2d] 
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AE  -  GC  »  4a!b2d2(e  -  X)  +  4a1a2cd(e  -  X)  +  8c(c  +  X)(e  -  X)(a1b2 

-  a2b i)  -  8bj^b2cd(c  +  X)  -  8a2bj^c2(c  +  X)  +  2(e  -  X) 

x  (d2  -  B2)(a2bj  -  a^b2)  +  2bjb2d(d2b2)  +  2a2bjc(d2  -  B2) 

-  4aia2cd(e  -  X)  -  4a^b2cd2  -  4aia2c2d  -  4a2bjd2(e  -  X) 

+  8bib2cd(c  +  X)  -  2b ^b2d(d 2  -  B2)  +  4a2bjcd2  +  8a^b2c2 
x  (c  +  X)  -  2a^b2c(d2  -  B2)  +  4a^a2c2d 

■  4d2(e  -  X)(a^b2  -  a2bj)  +  8c(c  +  X)(e  -  X)(a^b2  -  a2bj) 

+  8c2(c  +  X)(aj^b2  —  a2b^)  +  2(e  -  X)(d2  -  B2)(a2bj^  -  a^b2) 
+  2c(d 2  —  B2)(a2bi  -  a^b2)  +  4cd2(a2bj  —  ajb2) 

*  (a ib2  -  a2b i) [4d 2(e  -  X)  +  8c(c  +  X)(e  -  X)  +  8c2 

X  (c  +  X)  -  2(e  -  X)(d 2  -  B2)  -  2c(d2  -  B2)  -  4cd2] 

-  (aib2  -  a2bi){4d2(e  -  X  -  c)  +  8c(c  +  X)(e  -  X  +  c) 

-  2(e  -  X)(d2  -  B2)  -  2c(d2  -  B2)  -  4cd2] 

*  (aib2  -  a2bi)[4d2(e  -  X  -  c)  +  8c(c  +  X)(e  -  X  +  c) 

-  2(d2  -  B2)(e  -  X  +  c) 1 

m  (alb2  “  a2blM(®  -  X  +  c)(8c2  +  8cX  -  2d 2  +  2B2) 

+  4d2(e  -  X  +  c)  -  8d2c]  *  (ajb2  -  a2bj)[(e  -  X  +  c) 
x  (8c2  +  8cX  -  2d2  +  2B2  +  4d2)  -  8d2c] 

■  (a ib2  -  a2b i)  [ (e  -  X  +  c)(8c2  +8cX  +  2d2  +  2B2)  -  8d2c] 

-  2(a ib2  -  a2b1){(e  -  X  +  c)(4c2  +  4cX  +  d2  +  B2)  -  4d2c] 


Using 
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e  -  (1  -  r2)2  +  4r2  cos2  arc 


e  -  X 


+  c  =  (1  -  r2)2  +  4r2  cos2  a>x  -  (1  +  r2)2  +  4r2  cos2  an 


2  2 
+  4r  cos  arc  +  4r‘ 


1  -  2r2  +  r4  +  8r2  cos2  arc  -  1  -  2r2  -  r4  +  4r2 


2  2 
8r  cos  arc 


4c(c  +  X)  ■  16r2[4r2  +  (1  +  r2)2  -  4r2  cos2  an] 

=  16r2[4r2  sin2  an  +  (1  +  r2)2] 
d2  +  B2  *  16r2(l  +  r2)2  cos  an  +  16r2(  1  -  r2)2  sin2  an 
=■  16r 2(  1  +  r2)2  cos2  an  +  16r2(l  +  r2)2  sin2  an 


4  2 

64r  sin  an 


16r2(l  +  r2)2  -  64r4  sin2  an 


4c(c  +  X)  +  d2  +  B2  *  16r2[4r2  sin2  an] 


32r2(l  +  r2) 


(e  -  X  +  c)(4c2  +  4cX  +  d2  +  B2)  -  4D2c 

-  8r2  cos 2  an(32r2)(l  +  r2)  -  4 [ 16r 2( 1  +  r2)2 
x  cos2  an](4r2) 


256r4(l  +  r2)2  cos2  an  -  256r4(l  +  r2)2  cos2  ojx 


2\  2 _ 2 


Therefore, 


AE  -  GC 


(D-l) 


r  AE 
C  G 
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